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OS ! abstract 
OS 

Cosmological solutions of the Brans-Dicke theory are investigated by including a 
qj \ quantum effect coming from 1-loop correction of matter fields that couple to the 

scalar field. As the most serious result we face a cosmological "constant" in the 
original conformal frame which is shown to be "physical" after a careful analysis 
of the time variable employed in any of the conventional approaches. We find 
an "attractor" solution featuring no expansion during the radiation-dominated 
eras. To evade this unacceptable consequence, we suggest to modify one of the 
fundamental premises of the model, rendering the scalar field almost "invisible." 
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Among many versions of the scalar-tensor theory of gravitation, the prototype 
Brans-Dicke (BD) model JIJ is unique for the following four assumptions on the 
scalar field: (i) a nonminimal coupling of the simplest form; (ii) masslessness; 

(iii) no self-interaction; (iv) no direct matter coupling. 
Although the model may not be fully realistic, it still seems to deserve further 

scrutiny as a testing ground of many aspects of wider class of the scalar-tensor 
theories. The model has been studied, however, mainly as a classical theory. 
We attempt to take quantum effects due to the matter coupling into account. It 
has been argued that composition-independence that entails from the assumption 

(iv) above would be violated as a quantum correction [^]. We came to realize, 
however, that the suspected contribution is canceled by other terms arising from 
regularization;^ WEP is in fact a well-protected and robust property of the model 
beyond the classical level. 

We discuss in this note another quantum effect which, as it turns out, has 
serious cosmological consequences, beyond the extent of remedy expected by ad- 
justing the fundamental parameter of the model. For a possible way out we 
suggest to modify the assumption (iv), making the scalar field almost "invisible," 
thus reconciling with the absence of experimental evidences, still playing a cos- 
mological role. We also discuss how the scalar field acquires a nonzero self-mass 
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2 More details on this analysis will be reported in another publication. 
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due to the matter coupling, even having started with the classical assumption 
(ii), but likely in an entirely insignificant manner in practice. 
We start with the basic Lagrangian 



C x /^7 ( l -F{4>)R - e-g^d^d^ + L m ) . ( 1 ! 



with 

F{4>) = £0 2 - (2) 

We use the unit system of 8nG = 1.0 The scalar field and the constant £ are 
related to the original notation (p and lu, respectively, by 

<f = ^0 2 , and u = — . (3) 

We also allow e = ±1, a minimum extension of the original model to avoid an 
immediate failure. As we see shortly, e = — 1 does not necessarily imply a ghost 
in the final result. 

As the matter Lagrangian we choose, according to the assumption (iv), 

L m = -^(p+m )ij, (4) 

where ip stands for a simplified representative of the (spinor) matter field. 

In spite of (iv) couples in effect to the matter field in the field equation. 
This is inconvenient, however, when we try to apply the conventional technique 
of quantum field theory to the -matter coupling. For this reason we apply a 
conformal transformation such that the nonminimal coupling is eliminated ||: 

1 /2 

= fT 2 # w , with Q = (^0 2 ) . (5) 
Notice that we should have 

f > 0, (6) 

in order to ensure that the sign of the line element remains unchanged. 
We in fact find that (ID) is now put into the form 



£ = [z R * - <V<V + L, m J , (7) 

expressed in terms of the new metric g*^ and the field "0* = fi~ 3 / 2 -0. Also the 
canonical scalar field is now a as defined by 

0(a) = 0o eP°, (8) 



3 The units of length, time and energy are 8.07 x 10 _33 cm, 2.71 x 10 -43 sec, and 2.43 x 10 18 
GeV, respectively. Notice also that the present age of the Universe is ~ 10 60 . 
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where 

0- , 1 . - { ^~ (9) 
We emphasize that a is not a ghost if 

P 2 > 0, (10) 

even if e = — 1 Q; "mixing" between the scalar field and (spinless part of) the 
metric field provides sufficient amount of positive contribution to overcome the 
negative kinetic part. The condition fllCf ) is equivalent to 

3 

e^T 1 > -6, or euj > --. (11) 

We now have a direct a-matter coupling expressed in terms of the interaction 
Hamiltonian to which usual perturbation method is readily applied. We point 
out, however, that EEP, hence WEP as well, remains intact because the deviation 
from geodesic arising from the transformation @ is given entirely in terms of o 
which is independent of any specific properties of individual particles; the fact 
that any motion is independent of the mass, for example, verified in one conformal 
frame obviously survives conformal transformations. 

We call the conformal frames before and after the conformal transformation 
J frame and E frame, respectively.^ 



2 Quantum effect 

In E frame, in which we hereafter suppress the symbol * for simplicity, we consider 
one-loop diagrams as shown in Fig. 1, due to the interaction 

^ = mor 1/2 f 1 W#, (12) 

coming originally from the mass term in J frame. 

As will be shown, o may keep moving with time, and so does the mass. But a 
moves so slowly compared with any of the microscopic time-scales, that the mass 
m(t) at each epoch will be defined by 

m(t)=m C 1/2 r\o-(t)). (13) 

Now consider the 1 diagram (a). Its contribution is given by the potential 
of cr; 

V x {a) = im{a) f d A P Tr( 1 ) . (14) 
J \m(a)+ipj 

4 The name J frame is used, following Cho's suggestion ||, after P. Jordan B who was the 
first to discuss the nonminimal coupling. On the other hand, E frame is a reminder that this 
is a frame in which the standard theory of Einstein is formulated. 
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Our consideration will be restricted to the Universe which is sufficiently late to 
justify to ignore the effect of temperature and spacetime curvature. 5 

The integral in fll"4]) is quadratically divergent, but is expected to vanish if 
there is supersymmetry because the fermionic contribution given by (Jjjj) is can- 
celed by the same contribution from the bosonic partner. The cancellation would 
not be complete, however, if supersymmetry is broken at the mass scale 

M ssh = rm, (15) 

where the ratio r of M ss b to m, a representative of ordinary particles taken roughly 
of the order of GeV, would be ~ 10 3 -10 4 , which we naturally choose to be a true 
constant. 

The result may be given by 

V x {a) = C ir 2 m 4 = y 1 (0)e- 4/3,J , (16) 

where 

\/ 1 (O) = C 1 r 2 (m o r V2 0o 1 ) , (17) 

with C\ most likely of the order 1. Its sign, however, may not be known precisely 
because it depends on the details of supersymmetry breaking. If C\ > 0, ( pf) 
gives a positive exponential potential that would drive a toward infinity. We 
assume this to be the case. 



3 Cosmology 



3H 2 = 1 -^ + V + p, 

a + 3H& + V'{a) = 0, 
p + AHp = 0, 



(18) 

(19) 
(20) 



We now consider the cosmological equations with a classical potential V = V\ 
which is the only potential due to the assumption (iii): 



where p is the matter density which we assume, for the moment, to be relativistic. 
We also ignored possible terms representing the coupling between a and p. This 
would be justified for our purposes as long as we consider late epochs during 
which the coupling is sufficiently weak 0]. 

With V — Vx, a set of analytic solutions of flTB|)-(PUD are obtained: 

a(t) = t 1/2 , (21) 

5 The temperature will be lower than ~TeV, for example, if t>,10 -13 sec, much earlier than 
the epoch of nucleosynthesis. Spacetime curvature will be important only for £<,lCP 28 sec, 
corresponding to the temperature of 10 11 GeV. 
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a(t) = -^lnt, (22) 



2/5 

p{t) = \{l-\^)r\ (23) 



with 



16/3 

Notice that the condition p > is met if 

/T 2 < 4, or e£ _1 < -2, (25) 

which, combined with (^|), is satisfied only if e = — 1. This is the reason why we 
decided to allow an "apparent" ghost in J frame. Then fl2"5]) translates into 



r 1 >2, or cu>i (26) 

a much milder constraint than those derived from the observation. With e = — 1, 
however, ([TT| ) implies 

r X <6, or u<-, (27) 

which is ruled out immediately by the solar-system experiments, giving u;,> 10 3 
0. We nevertheless continue our analysis as long as theoretical consistency is 
maintained. 

From also follows 

<P(t) = U Xj \ (28) 

where 0o is determined by identifying ( |TTD with (pi]); 4 = 16Cir 2 m 4 /3 2 £~ 2 , which, 
if used in flilf ), yields 

m(t) = i(C 1 r 2 /3 2 )" 1/4 r^. (29) 
It is interesting to notice that the behavior m(t) ~ t -1 / 2 follows simply because 



the potential should be proportional to m 4 , as shown in fllBT ), and it must decay 
like t~ 2 because it is part of the energy density appearing on the right-hand side 
of the 00-component of the Einstein equation. 

Obviously the assumption r = const is crucial in the above argument. We 
could obtain m = const if r(t) ~ but with a highly unreasonable consequence 
that r should be as large as ~ 10 63 at the Planck time. Also the dependence 
V ~ m 4 is common to any diagrams of many _1 's, as in (b) and (c) in Fig. 1. 
Including them results simply in affecting the overall size of the potential. 

We arrived at ( p9|) in E frame in which the standard technique of quantum 
field theory can be applied. We should also notice, however, that we use some 
type of microscopic clocks in most of the measurements. The time unit of atomic 
clocks, for example, is provided typically by the frequency ~ ma 4 . It is also 
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important to recognize that the cosmic time is usually assumed to be measured 
in the same time unit. If the time unit r(t) itself changes with time, the new 
time t measured in this unit would be defined by 

dt , „ 

dt = — 30 
r{t) 

In conformity with special relativity, the scale factor a(t) in Robert son- Walker 
cosmology is transformed in the same manner: 

a(t) 

T{t) 

These two relations can be combined to a conformal transformation 

ds 2 = r~ 2 ds 2 , or g^ u = r 2 g^ v . (32) 

In the prototype BD model, there is no mechanism to make a time-dependent, 
hence r ~ m _1 . Combining this with m ~ as derived from flilf ), and also 
comparing ( p2|) with (|5]), we find that (|32}) implies going back to the original 
J frame, as it should because it is the frame in which mass mo is taken to be 
constant. 

We now try to solve the cosmological equations in J frame, in which we also 
suppress tildes to simplify the notation. It is also interesting to find that the 
behavior V\ ~ 0~ 4 as indicated in ([16]) shows that this potential in E frame can 
be derived from a cosmological constant in J frame as given by 

A = C x r 2 ml. (33) 

The quantum effect computed in E frame amounts to introducing A back in the 
original conformal frame. The field equations in J frame are given by 

2i P G tlv = r^ + T / t-^A-2(^n-V M V,)^, (34) 
U<P = /3 2 (T-4A), (35) 
V M T^ = 0. (36) 

Notice also that is the matter energy-momentum tensor while 

T% = e - l -g, v {d4>?) • (37) 

Assuming spatially uniform 0, we derive the cosmological equations: 

6^H 2 = 6^<p 2 + A + p-6Hp, (38) 

<p + 3H<p = 4/5 2 A, (39) 
p + 4Hp = 0, (40) 
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where we have chosen T = confining ourselves to the radiation-dominated era 
even if A comes from nonzero m . 

As a heuristic approach, let us choose 

H = 0. (41) 

Then (SO) leads to 



p = const. (42) 
Using fl41|) in (^) we obtain Cp = 4(3 2 A, which allows a solution 

<p(t) = 2f3 2 At 2 + tprf + <p - (43) 

Notice that we have chosen A > hence C\ > so that a falls off the potential 
slope toward infinity in E frame. This implies that ip increases also in J frame 
as indicated in ([S]) if f3 2 > 0. This is the very condition, however, which is in 
contradiction with the observation, as already discussed in connection with (|27|). 
Taking aside this drawback for the moment again, we expect 

(44) 




at sufficiently late times. Using this together with ([II]) and ( f4"2|) in (|38| ) gives 

P = 3A^—^, (45) 
H 6£-l' v ' 

which is positive if (^) and Q2"7D are obeyed. 

An example is shown in Fig. 2, in which we see how H approaches zero, 
much faster than t~ l ; the Universe quickly becomes stationary after alternate 
occurrences of expansion (H > 0) and contraction (H < 0). This together with 
other similar examples indicate strongly that there is an "attractor" to which 
solutions of different initial conditions would approach asymptotically. Fig. 3, 
which is a 2-dimensional cross section of the 3-dimensional phase space of ip, ip 
and p, illustrates how different solutions are attracted to a common destination 
given by (|44|) and (f|^), which represent in fact a curve in the whole phase space 
as one finds because of the relation <p 2 = up 2 j (2^<p). A trajectory for a set of initial 
values proceeds along, spiraling around and coming ever closer to this curve. 

One might be tempted to compare our solutions with those in Einstein's model 
with a negative cosmological constant A < 0, but of course without 0. This 
model allows a static solution p = —A and H = 0, but the Universe would 
never become stationary in contrast to our solutions; p with a sufficiently large 
initial value decreases toward the minimum —A but bounces back to increase in 
a touch-and-go fashion. 
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In this way we come to conclude that the BD model corrected for an important 
quantum effect should result in a steady state Universe, which is totally unac- 
ceptable in view of the success of the standard model in understanding primordial 
nucleosynthesis. 

The constant scale factor in J frame may be interpreted also from the analysis 
in E frame, in which the length unit is provided by r ~ <\> ~ t 1 / 2 which increases 
in the same rate as a ~ t 1 / 2 shown in (plD |§ . We also find that (f> ~ t as given by 

1 /2 

(fPD in J frame and <p ~ t* in (P8"D, in which the symbol * is restored in E frame, 
are consistent with each other, since t ~ tV 2 is a consequence of the relation 
(l30|).p| These observations seem to support (|20|) which is only approximate unlike 
its counter part fl4"0"|) in J frame. 

On the other hand, one may ask if there is any sensible solution with a ~ 

t l/2 

in J frame. In fl40|) we substitute 

H = (l/2)t~\ (46) 



thus obtaining 
Then (R9[) becomes 



P = Pot' 2 - (47) 



(p + = 4/5 2 A, (48) 



2 



which is solved asymptotically: 

4 




cp^jf3 2 At 2 , or 0^ t /e^(l-6/3 2 )At. (49) 



We then find that the right-hand side of (|38| ) is given by 

p + jj (3 - 16/3 2 ) A. (50) 

Now from (|46D and (|49|), the left-hand side of (^) should be time- independent. 
This can be matched with the situation in which p given by (|47| ) decreases rapidly 
to be negligible compared with the second term of (|50| ) ; implying that the Universe 
becomes asymptotically "vacuum dominated," again an unrealistic conclusion. 
Even worse, ignoring p in the right-hand side of ([5(]) and using ( fi9"| ) on the left- 
hand side of (|38| ) yields 

P 2 = ~, (51) 

which on substituting into (H) gives 

r 1 = 0, (52) 

hardly a realistic result .[] 

6 Apply the replacement, t — * t, t — > t*. For the dust- dominated Universe with a* ~ i*^ 3 , we 
find a - t 1 / 6 . 

7 The same analysis applied to the dust matter results in e = —1 and £ = —2/3, being 
inconsistent with (|6|). 
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4 Discussions 



We add that our argument of choosing J frame is independent of whether we 
literally use atomic clocks to measure something during the epoch in question. 
It is simply in accordance with realistic situations that analyses are based on 
quantum mechanics in which mass of every particle is taken to be truly constant. 

We admit that there should be some other quantum effects to be included. 
The result obtained here is, however, so remote from what would be expected 
from the standard scenario that it is highly unlikely that those "other" effects 
conspire miraculously to restore the success in the nucleosynthesis, among other 
things. It seems that we need some more fundamental modification of the model. 

A possible way out is to abandon the assumption (iv) about the absence of 
the direct <fi coupling to the matter in J frame. As an extreme counter example, 
we may replace the matter Lagrangian @ in J frame by 

L^ = -$(#+/0)V, (53) 

where / is a dimensionless Yukawa coupling constant, ip is massless in J frame, 
while in E frame we obtain the mass m = f^," 1 ^ 2 (in units of (87rG) -1 / 2 ) which 
is independent of a. The scalar field is decoupled from if) in E frame, hence is 
left invisible through the matter coupling; it plays a role only in cosmology, most 
likely as a form of dark matter. With constant mass in E frame, which is now 
physical, we may reasonably adjust parameters such that the standard scenario 
is reproduced to a good approximation. 

There might be intermediate choices between the prototype model and this 
extreme model. Then we may expect the matter coupling generically much weaker 
than 

H' ai = -Pm(tfH><r, (54) 

in the prototype model, hence evading immediate conflicts with the test of WEP 
and the constraint from the solar-system experiments.^ Needless to say, G is 
predicted to be constant, by construction. 

With modifications of this nature in mind, we add a comment on the mass 
term of the scalar field, which is ought to arise from V\ as given by fll6P . We 
should be interested here in a fluctuating component o~\ (x) which is responsible 
for the force between local mass distributions, to be separated from the spatially 
uniform component o~o(t) evolving as the cosmic time t; 

a(x) = cr (t) + <7i(x), (55) 

which satisfies 

na _ V \a) = 0, (56) 

from which (|l||) derives. 

8 See Ref. 0] for a model of this type. 
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The cosmological component do is a solution of ( |56"D with a\ dropped; 



□ (7o + 4^!(0)e- 4fe =0. (57) 
If we use this in (|56"D for the entire field (|55"|) , we obtain 

□ a! - 4/3Vi (0)e- 4/3<TO (l - e~ 4/3<Tl ) = 0. (58) 
Expanding the terms in the last parenthesis, we find 

nai = fiVi + • • • , (59) 



where, using (|22| ) for cr (t)|] 

/i 2 = 16/3V!(0)r 2 , (60) 

which is V"(cr) at o = cr . 

This shows that the scalar field does acquire a "mass" even though the poten- 
tial has no stationary point, but the range of the force mediated by o~\ is basically 
given by t, which is the size of the visible part of the Universe at each epoch. The 
force-range at the present epoch can be as "short" as 10 5 ly if 16/3 2 Vi(0) ~ 10 10 , 
in contrast to (p4|). It is rather likely that the force can be considered to be 
infinite-range in any practical use. 

We may relax the assumption F ~ cf 2 in the prototype model. We recognize, 
however, that the relation of the type V(a) ~ m 4 (<r) as in (|T6| ) is quite generic 
and so is m ~ t -1 / 2 according to the argument following (^). This makes the 
conclusion (^TJ) almost inevitable, as long as the assumption (iv) is maintained. 

As another aspect of more general F((f)), we point out that the factor £ _1 / 2 -1 
in (TT) is in fact F~ x l 2 . It then follows that the potential as given by (|Ttj|) 



generalizes to 

V\{p) ~ m 4 ~ F~ 2 . (61) 
The relation (M) is also traced back to 



We may then expect that the potential as a function of o would have a minimum 
if F((p) has a maximum, the same conclusion as in Refs. 0. 

The potential minimum should act, however, as an effective cosmological con- 
stant, which might present another serious conflict with observations unless it 
remains below the level of ~ to 2 ~ 10 -120 at the present epoch. In this respect 
we have an advantage in the potential having no stationary point. 

As the last comment we point out that the occurrence of a ghost which was 
required to give positive matter density is not entirely unnatural from the point 

9 For dust-dominated Universe, the right-hand side is doubled. 
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of view of unified theories. Consider, for example, Kaluza-Klein approach to 4 + n 
dimensional spacetime. The size of compactified n-dimensional space behaves as 
a 4-dimensional scalar field, which is shown to have wrong sign in the kinetic 
term. This model provides also one of the natural origins of the nonminimal 
coupling. 
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(a) (b) (c) 

Figure 1: Some of the one-^-loop diagrams, shown only up to the third order in 
Each blob represents -1 appearing in (12), a collection of many vertices of 
cr's, while a solid line is for the ip line. 
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Figure 2: A solution of (38)-(40) with the initial values ip = 0.25. (p = 0.0, 
0.1 at hit = 1. We chose A = 0.5 and f = 0.4. 
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Figure 3: Trajectories (lnt = 1 — 30) of the solutions of different initial values y?o 
and po) as shown in the parentheses, with other parameters the same as in Fig. 
2. The point of convergence (0 2 = 1.42857, p = 0.21429) is given by (44) and 
(45). The dotted line is for H = 0; its left-side for H > 0. 
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